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We construct a generic coarse-grained dynamics of a thin inflexible planar layer of polar-ordered 
suspension of active particles, that is frictionally coupled to an embedding isotropic passive fluid 
medium with a friction coefficient F. Being controlled by F, our model provides a unified framework 
to describe the long wavelength behaviour of a variety of thin polar-ordered systems, ranging from 
wet to dry active matters and free standing active films. Investigations of the linear instabilities 
around a chosen orientationally ordered uniform reference state reveal generic moving and static 
instabilities in the system, that can depend sensitively on F. Based on our results, we discuss 
estimation of bounds on F in experimentally accessible systems. 


I. INTRODUCTION 

The emergence of large-scale collective dynamics is one of the most intriguing and fascinating features of a large 
variety of driven, active systems made of active particles [l|. These are generally elongated and their direction of 
self-propulsion is set by their own anisotropy (i.e., the two ends are distinguishable, hence polar), instead of being 
determined by an externally imposed field. In contrast, active nematics 0, made of active particles which are head- 
tail symmetric, do not show any self-propulsion. These active systems, polar or nematic, are generically characterised 
by the existence of orientationally ordered states. These are nonequilibrium analogues of the equilibrium nematic 
liquid crystals. There are numerous examples, which include both living systems {living matter) as well as their 
artificially prepared non-living analogues. Biological examples of active systems include both small and large in-vitro 
and in-vivo systems, e.g., reconstituted bio-filaments and the associated motor proteins the cytoskeleton of living 
cells and bacterial suspensions 0, cell layers and also larger-size objects, e.g., flock of birds or school of fishes 0. 
Analogous non-living examples of active matter systems also arise in various contexts, e.g., layers of vibrated granular 
rods (3 and colloidal or nanoscale particles propelled through a fluid by catalytic activity at their surface All 
these examples of active systems are distinguished by a local energy supply in the bulk that drives the systems away 
from equilibrium. This is in contrast to other well-known examples of driven systems, e.g., sheared systems, where 
the external drives act at the boundaries. For instance, in cell biology contexts, this supply of energy takes place due 
to the hydrolysis of adenosine triphosphate (ATP) to adenosine diphosphate (ADP) and other phosphates (Ph) by 
the molecular motors, thus converting the chemical fuel into mechanical motion. 

Despite hugely varying microscopic details, different active matter systems display a host of intriguing nonequilib¬ 
rium phenomena with generic features independent of system details, e.g., pattern formations, wave propagations, 
oscillations and unusually strong fluctuations iEl. Due to the large number of diverse microscopic variables present 
(especially in the cell biology context), the level of complexities in active matter systems at microscopic levels is very 
high. Instead, it is convenient to formulate the coarse-grained dynamics of the active systems based on identifying 
global features, e.g., the presence or absence of conservation laws, symmetries, the presence of appropriate broken 
symmetry variables and the nature of the underlying momentum damping. These are similar in spirit and nonequi¬ 
librium generalisation of the general principles and laws developed to describe the statistical mechanics and dynamics 
of the ordered phases in equilibrium systems [l^. These active fluid theories, parametrised by a set of phenomeno¬ 
logical constants ilIM3 , serve as as generic coarse-grained descriptions for a driven orientable fluid with nematic 
or polar symmetries and are particularly useful to uncover and elucidate the long wavelength behaviour observed in 
very different physical systems and at very different length scales iiEIIl. 

In a bulk fluid (both active and passive) the viscosity damps out any local momentum gradient and thus reduces 
any relative velocities between neighbouring regions. The total momentum of the system is however kept conserved; 
such systems are known as wet active matters in the language of Ref. [l^ : see, e.g. Refs. @[11. In contrast, for 
systems resting on a rigid substrate (e.g., a layer of active fluid on a solid substrate) there is a drag on the system 
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acted typically through a no-slip boundary condition on the active matter velocity at the active matter-rigid substrate 
interface. This drag leads to nonconservation of the momentum of the active system and cuts off any long-ranged 
hydrodynamic interactions. These are known as dry active matter in the classification used in Ref. M and have been 
studied extensively, see, e.g., [i,a[ia,i3- The properties of active matter systems are often considered in the form 
of thin, quasi two-dimensional (2D) layer. Such quasi-2D active matter systems exist both in-vivo and in-vitro: cell 
cortex [m or the cortical actin layers and cell ruffles, e.g., lamellipodia are examples belonging to the former 
category, where as reconstituted actin layers on liposomes [23l| are examples of 2D in-vitro active fluid systems. 

Inspired by the current studies on both wet and dry active matters and their significant differences in terms of 
their long wavelength properties, we study a generic 2D polar active matter layer, where the active particle system 
is embedded inside a three-dimensional (3D) bulk isotropic passive fluid. The active fluid and the embedding passive 
fluid interact via a mutual friction at the interfaces of the active fluid-bulk fluid interface, leading to momentum 
damping of the active particles. To this end, we construct a set of 2D continuum equations of motion for the local 
orientation and number density of the polar active species. Our model, parametrised by the interfacial friction T, 
provides a unified coarse-grained description of the dynamics of polar ordered 2D wet and dry active matters and free 
standing 2D films. In a linearised treatment about a chosen orientationally ordered uniform state, we find the linear 
instabilities in the system. We also study the nematic limit of the dynamics. The nature of the linear instabilities 
are found to depend sensitively on the magnitude of F relative to the viscous damping. Our results may be used to 
estimate bounds on T in possible physical realisations of our model, e.g., reconstituted actin filaments deposited on 
a liposome embedded in a fluid medium. In addition, in an in-vivo system of two eukaryotic cells with a substantial 
area of contact, the dynamics of the cortical actin layers of the two cells on both sides of the contact plane should be 
describable by our dynamical equations at a coarse-grained level. Nonetheless, our formulation is sufficiently general 
and does not specifically relate to any particular cell biological example. The rest of the article is organised as follows: 
In Sec. m we define our model and set up the basic equations of motion. Then in Secs. lIIlXllIIIBI and llllCl we analyse 
the instabilities for high, intermediate and low values of the mutual friction. Then in Sec. II VI we briefly compare the 
linear instabilities in the different regimes of the model, delineated by the magnitude of the mutual friction. In the 
next Sec. |Vl we analyse the nematic limit of our model dynamics. We discuss and summarise in Sec. IVII Finally, we 
provide some calculational details and then obtain the ambient velocity profiles in the Appendices. 


II. MODEL EQUATIONS 

We consider an inflexible thin planar layer of a viscous active fluid with a vanishingly small thickness, located at the 
xy-plane, i.e., at ^ = 0. We treat it as a quasi 2D system, for which a 2D description should be appropriate. The local 
number densities of the active species and the solvent are p(x) and ^i(x), x = {x,y), respectively. The active fluid 
layer, with a 2D viscosity 77 , is embedded in a 3D passive incompressible ambient fluid with a 3D viscosity 77 ', both 
above {z > 0) and below (z < 0); see Fig. [T] for a schematic diagram of our model system. It is not unusual to treat 
thin active fluid layers as quasi-2D systems; see, e.g.. Ref. [ 2 ^. We expect this 2D description with a 2D viscosity to 
be good for really very thin system such that any variation of the physical quantities along the thin direction may be 
neglected. 

The centre of mass velocity of the active particles and the solvent combined is given by v. The total number of 
both the active and solvent particles are separately conserved: The continuity equations for p and (/) are written as 

dtp + V-Jp = 0 , ( 1 ) 

dtc^ + V-J^, = 0 . ( 2 ) 

Here, V = xdjdx -I- yd/dy is the 2D gradient operator, x, y are the unit vectors along x- and ^-directions. The 
particle currents Jp and can be expressed in terms of the 2D centre-of-mass velocity v and the diffusion current j. 

Jp = pv-bj, (3) 

J0 = (Av-j. (4) 

Here, the molecular masses of both the active and solvent particles are assumed to be equal to unity for calculational 
convenience. We are interested in an orientationally ordered state of the model system. To this end, we introduce a 
2D local polarisation vector p = {Px,Py), with a fixed magnitude, p^ = I, as appropriate for an orientationally ordered 
state. Microscopically, it describes the local orientations of the actin filaments or bacteria. We consider the active 
fluid to be overall incompressible, i.e., V • v = 0. Our chosen reference state is defined by Pa, = 1 with no macroscopic 
overall flow, i.e., {va) = 0. Note that this does not rule out finite velocity of propagation (or a nonzero macroscopic 
current) of the active particles (see below). In the Stokesian limit of the flow dynamics, the force balance equation 

VcCTa/g — dpii -b F /3 = 0, with a, (3 = x,y (5) 
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FIG. 1: (Color online) Schematic diagram of our model active fluid layer spread along the xy-plane. Arrows indicate polar 
active particles, aligned predominantly along the i-direction (see text). 


yields the generalised Stokes equation for v. The 2D pressure 11 may be eliminated by using the incompressibility 
condition (V ■ v = 0). Here, aap is the total stress tensor and external forces Fp are the tangential stresses of the 
embedding fluid on the two sides (top and bottom) of the active fluid layer 

Fp = ■q'idWp + ( 6 ) 

where e —>■ 0; v'(r) (with r = {x, y, z)) is the 3D ambient fluid velocity. 

ambient fluid 



active fluid layer 


FIG. 2: (Color online) Schematic diagram of our model active fluid layer spread along the a:j/-plane (« = 0). The interfacial 
friction acts at z = ±e, e —>■ 0 (see text). 

In the spirit of linear response theories [l^, the dynamics of the active fluid layer is described in terms of linear 
relations between the thermodynamic fluxes (cr®«, ja,Pa) and the corresponding generalised forces {vap,dajl,ha) [^, 
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[l3, m, [2^. Here, cr®^ is the symmetric part of the deviatoric stress 

</3 = + PtVaVp - alp, (7) 

with alp = {pahp—ppha)/2 is the antisymmetric part of the stress tensor, ha being the thermodynamic force conjugate 
to polarisation pa- Further, pt = p + 4> is the total density of the two species combined and Vap = [daVp + dpVa)/^ is 
local strain rate tensor. The term ptVaVp is the Reynold’s stress in the active fluid. In addition, P is the convected 
co-rotational derivative of the polarisation vector given by 

Pa = = 9tPa + VpdpPa + UJaPPp, ( 8 ) 

with LOap = ^{daVp — dpVa) is the vorticity tensor. Furthermore, p, = Pp — pp is the effective chemical potential; 
Pp and pp are individual chemical potentials of the active particles and the solvent molecules, respectively. For 
simplicity we consider the dilute limit of the active particles p (f) ov pt ~ (j), i.e., ~ (pv. In this limit, the overall 

incompressibility (which implies pt = const.) is equivalent to considering (j) = const., so that the dynamics of (j) can be 
neglected and we consider the dynamics of p alone. In this dilute limit, p may be replaced by the chemical potential 
Pp for the active particles. 

The stress field is assumed to contain a nonequilibrium active stress of the form 

CTaf = C{p)^PPaPp. (9) 

Microscopically, cr“^‘ is due to the local nonequilibrium dynamics of the active particles. The coarse-grained form 
(0 may be obtained Q by noting that the force applied by an active particle on the fluid surrounding it is same as 
that applied by the fluid on it, considering the total forces exerted by a collection of active particles with a given 
centre, each exerting a point force proportional to and parallel to ±p and expanding it up to the lowest order in 
spatial gradients. Therefore, the magnitude of should depend on the local density p of the active particles; hence 
the form Q' = C,'{p)- We write p = po + ^Pi where po is the mean active particle density and 5p are fluctuations 
(assumed small) about pg. Expanding for small dp, we write C,'{p) = C + C<5p, where C = C^(Po) and Q = i9CV'5p|p=po' 
Parameter Ap represents the strength of a^p-, the latter is said to be contractile or extensile depending on whether 
Ap is negative or positive, respectively; Ap is a measure of the rate of supply of (free) energy that pushes the system 
out of equilibrium; in the context of the cortical actins in a cell, it is the hydrolysis of the ATP molecules to ADP and 
phosphates that supplies this energy; Ap = patp — Padp — pPh where patp, Padp, pph are the chemical potentials of 
ATP, ADP and phosphate molecules. Parameter Ap has the dimension of energy/(mass.mole). Numerical estimation 
of Ap is not easy: In the particular context of cell biology, one may use the fact that approximately 7 kCal energy 
released per mole of ATP due to its hydrolysis. Since 1 molar mass of ATP ^ 500, we obtain from its definition 
Ap ~ 7kCal/{500gm/10'^^), the free ener gy r elease per unit mass per molecule. 

The relevant linear flux-force relations |16l . , that include the active stress contribution to the stress and allow 

for polar terms, i.e., not invariant under p -A —p, are 

^aP = ‘^VVap+C'{p)ApPaPp + YiPahp+Ppha)-^{PadpPp+ppdaPp), 

ja = -'lppdahp + ^ha + KpPaAp + wApdp{pPaPp)-^PpidaVp+dpVa), 

Pa = — + XlPaAp + VippVaP “ \dajlp + A2 (p • V)paAp -|- XpdapAfJ.. 

7o 

Coupling constant i/i denotes the equilibrium flow-orientation coupling [^ : similarly, eg denotes symmetry-allowed 
equilibrium couplings between the flow and the particle current [25| . Parameter jpp > 0 is a mobility coefficient 
(an equilibrium coupling constant) and related to the diffusion coefficient. In addition, particle current jp should 
have active contributions /CpApp and wApdp{ppaPp), such that there should be an active macroscopic current of the 
particles in the direction of p, with amplitudes proportional to Ap. In addition, A is a cross-coupling equilibrium 
coupling constant. (In general, A may be a tensor reflecting the anisotropy of the polar ordered state; we neglect this 
here.) Notice that in (|12[) we include two symmetry-permitted active terms with coefficients A 2 and Ap, respectively; 
the A 2 -term is a self-advection term, (not considered in Refs. [H, [ 2 ^). Since the active particles tend to display 
macroscopic motion with respect to the embedding fluid even in their fully ordered state (no distortion), microscop¬ 
ically the A 2 -term represents advection of the local distortions in p by p. The Ap-term is a nonequilibrium partial 
pressure term, modelling motion of the active particles along or opposite to the concentration gradient (depending 
upon the sign). Coefficients , Kp, Xi, X2,w and Ap are ’’active coefficients”, i.e., coefficients of different active terms 


( 10 ) 

( 11 ) 
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in Eqs. (I10m2l) . Out of these, Kp,X 2 and Xp are coefficients of the different polar terms, which break the symmetry 
under p —p, where as w, ^ and ^ are coefficients of the nematic active terms in the dynamics. Thus, in the nematic 
limit of the model, Kp, A 2 and Xp are all zero, and the only source of nonequilibrium drive is the active stress ([9]) and 
the active particle current represented by the w-term. For reasons similar to the p-dependence of active coefficient 
Kp should depend on p. We write for small density fluctuations @ 

d K 

^p(p) = ^0 T ^pp — clp lp=P0 5 (13) 


where kq = Kp(po) depends on the mean density and Kpp incorporates the effects of the fluctuations of p about po- We 
ignore any p-dependence of Xp and A 2 and treat them as a constant, since we are interested in a linearised treatment. 
Thermodynamic forces h and Pp are defined as follows 


h 


Ol 


5 ^ _ ^ 

6pa ’ 5p 


(14) 


where is a free energy functional that controls the relaxation of the system to its thermal equilibrium state in the 
absence of any activity. At the bilinear order in fields 


^0 = 1 + A{5pf + 2xpV • p]. 


(15) 


where H is a 2D Frank elastic constant (we have assumed equal Frank’s constants for simplicity), A ^ Tpo is an 
osmotic modulus with T being the temperature when the system is in thermal equilibrium, x provides a symmetry- 
allowed coupling between the density fluctuations and splay. Assuming the minimum free energy configuration to be 
given by a uniform configuration p = const, and p = pa everywhere, we must have < AD. From (1151) . we find 
ha = = DX/'^pa + X^aP and Pp = = A5p + xV • p. Eliminating 11 and using the forms of h and Pp, we 

obtain to the lowest order in spatial gradients (see Appendix |A| 

-h (ApPp^d^p^ + (ApPp^dyPy + (ApPp^d^p = -Fp, (16) 


where we have linearised about Px = p = Po. 

For an isotropic, passive ambient fluid medium, in the low Reynolds number limit and for small masses, its velocity 
v^, i = x,y, z satisfies the Stokes Eq.: 

p'vlvl = V3.n', (17) 

valid for both the super- {z > 0) and sub- {z < 0) phases, 11' is the ambient fluid pressure and V3 is the 3D gradient 
operator. The boundary conditions on are as follows: 

• No flow at infinity: At both 2 : —>■ ± 00 , n' should vanish. 

• Balance of the normal stresses of the ambient fluid at the 2D active fluid layer, 

• Due to the assumed inflexibility of the active fluid layer, the normal velocity of the ambient fluid at the active 
fluid layer should be zero: v'^{z = e) = 0 = v'^{z = —e). 

• Boundary conditions at the active fluid-bulk fluid interfaces requires careful consideration; see Fig [5] The most 
common boundary condition used in this context is the ” no-slip” condition, i.e., equality of the active fluid 
velocity and the in-plane component of the 3D ambient fluid velocity at the top and bottom interfaces between 
the ambient fluid and the active fluid layer. We generalise this by allowing a slip. We implement this by 
introducing a slip coefficient of friction^ such that the shear stresses are balanced by the friction forces at the 
interfaces. This implies (using v'^ =Q z = ±e) 

= T{v'p\^-vp), (18) 

= -T{v'p\-^-vp)ffi = x,y (19) 

where T is the slip coefficient of friction at the upper and lower interfaces (we assume equal friction at the upper 
and lower interfaces for simplicity); this allows us to define a slip length Ig ^ pV^- 
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Notice that for a finite F, boundary conditions (l20l) implies partial slip between U/ 3 (a;, y) and v'p{x, y,z = ±e), f3 = x,y. 
While the no-slip boundary condition is more conventionally used, on mesoscopic scales, however, instances of violation 
of the no-slip boundary conditions are known. For instance. Ref. [2^ has shown that beyond a critical shear stress that 
depends strongly on the surface roughness, departure from the no-slip conditions may be observed. It has also been 
found that upon addition of surfactant in the fluid, the boundary condition changes from no-sfe to partial slip [2^ . 
In addition, there are now strong evidences in favour of slip in polymer melts; see, e.g. Refs. [30|, |3lJ- Furthermore, 
it has been demonstrated in Ref. how a large slip at a liquid-liquid interface may be introduced experimentally. 
Friction has been considered in various active fluid flow problems as well; see, e.g.. Refs, for various theoretical 

and expreimental studies. While no systematic measurements of slip at interfaces involving active fluids are known, 
the above existing results suggest that considering the complex internal structure of the active fluid (e.g., the presence 
of actin filaments), a partial slip at the interfaces between the active fluid layer and the 3D embedding fluid cannot be 
ruled out. Recently, it has been shown that a significant reduction of the sliding frictional forces between two bundled 
F-actine filaments may be achieved by coating the F-actins with polymeric brushes [s^. Thus it is important to study 
implications of finite slips in an active fluid problem in a simple set up, which we set out to do below by using our 
model system. Notice that vanishing ui, at z = ±e implies that the shear forces F/s on the active fluid layer as given 
in ([n|) take the simpler form 


Fp = r]'{dzVp\^=e - dzVp\z=-e). ( 20 ) 

By using the boundary conditions prescribed above, together with the incompressibility of the ambient fluid Va-v' = 0, 
Stokes’ Eq. (flTl) may be solved to yield i = x,y, z (see Appendix) and obtain F)?. 

Equation of motion for the orientational field Pa may be written combining equations ([HI) and ca). 

dtPa + VpdfjPa + UJapPp = — + XiPaAfJ. + VippVaP “ XdaPp -|- A2(p • V)paAp + XpdapAp. (21) 

7o 

With px = 1 dehning the reference state, py is a broken symmetry (slow) mode. We linearise (1211) above about px = 1 
for small py. This yields 

dtPy = X2ApdxPy + XpApdyp + -f ^ dxVy - XAdyp - Xxdlpy. (22) 

7o 2 2“ 

Note that in Eq. (I22[) . p enters into the dynamics of Py through both equilibrium and nonequilibrium contributions. 
Both are equally relevant being the lowest order terms in gradient expansions. 

The equation of motion for p is obtained by using Eqs. (P, m and (ED. Up to the order the equation of motion 
for p in the Fourier space, linearised about = 1, is obtained as (set A = 1) 

dtp = -Ippff'P + wApqlp -I- wAppoqxqyP - iApKoqyPy - iApUppqxP + iDXqyq^Py -|- (23) 

where q = (g^,, qy) is the in-plane Fourier wavevector, conjugate to x = [x, y). 

Notice that in the linear equations dT51) , (12^ and (E^ there are seven active coefficients (excluding Apt) , which are 
introduced in the standard active fluid models Out of these, A 2 ,Ap,Ko and Kpp control the conditions 

for instabilities (along with the sign of Apt) for both the high friction and intermediate friction cases (see below). 
In terms of an underlying equivalent agent-based microscopic dynamics, we expect all these coefficients to depend 
upon the local density of the active particles and the specific alignment rules (favouring nematic or polar alignment). 
Thus, it is reasonable to expect that all the seven active coefficients are not independent parameters. On dimensional 
ground we argue that the two nematic active coefficients C, and ^ should be related as C ~ poC the pairs of polar 
active coefficients in the active particle current (kq) ^pp) and in the active alignment (A 2 , Ap) are related as KppPo i^p 
and ApPo A 2 , respectively. With the expectation that the polar alignment polar and active current terms originate 
from same underlying (polar) microscopic rules, we expect them to be mutually simply related. Again on dimensional 
ground we expect A 2 Kpp. Note however that in the all the above heuristic relations, there are dimensionless 
proportionality constants which we cannot obtain on simple physical ground. We would like to emphasise that all 
these parameters are just phenomenological constants, similar to the parameters which appear in the continuum 
theories to describe the statistical mechanics and dynamics of the ordered phases in equilibrium systems [H and 
cannot be calculated within our theory. It should in-principle be possible to relate these coefficients to and calculate 
them from the specific microscopic rules for agent based models for active systems; see, e.g.. Refs. (H-Hlj. These are, 
however, outside the scope of the present study. In what follows below, we ignore this issue for simplicity and treat 
all the seven coefficients as independent model parameters. In addition to these seven active coefficients, the friction 
coefficient F is not an active coefficient. This enters into the dynamics through the boundary conditions and is a 
model parameter that has been introduced by us and is central to the present discussion. To our knowledge, no good 
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estimate about the magnitude of F is available; we therefore treat F as a free parameter in our model. In general, 
Eqs. (froi) . (1^ and (1^ can be solved in principle for arbitrary values of F. Nonetheless, it is instructive to consider 
three different limits of F, characterised by Isjr]', system size L and thickness d of the 2D active system, and analyse 
them separately, as discussed below. 


A. High friction limit 


We consider a ’’large” F: F v'/L- Formally, we consider the dynamics in the limit F —oo (equivalently, —>■ 0); 
this is valid for wavectors qlg -C 1, or for a system of linear size L, Ig/L 1. Thus, the system size must be much 
larger than the slip length. The stress balance equations dn and (jl9|) yield 

= v'aU^-e = Va, a = x,y. (24) 


Thus, there is no slip between the ambient fluid velocity at the active fluid layer v'J\z=±e and the active fluid velocity 
Va ■ Equation IHM forms one of the boundary conditions on the ambient fluid velocity v '^. 

Forces (shear stresses) F^, then may be expressed as (see Appendix [B) see also Ref. j4^), 


Fj; = -2qvx'n', (25) 

Fy = -2qvyr]', (26) 


where q = (q^^qy) is the in-plane Fourier wavevector. Putting the values of (l25ll and (l26)) in the Stokes equation (fTHI) . 
the expressions for Vx and Vy can be derived up to the lowest order in q linearising about px = 1 and p = pq. 


V 


V 


X 


y 


Mldy 

2r]'q^ 

* 2-q'q^ 


Cq^ 

^PPy + ^^^^PPy 


2ri'q^ 


App, 


Mldy 

2ri'q^ 


App. 


(27) 

(28) 


Thus, Va at 0{q^) has only active contributions. 

Equation (1221) may be written by substituting for Vx and Vy from Eqs. (I27t and (l28l) . We thus obtain 


dtPy = 


-Dq^Py + ixqyP 
70 


Ari'q 


[{i^i - l)qy - {vi + l)qj 


iXpApqyp + i\2ApqxPy - iXqyP - Xxq^Py 

(Apqxqy 




Py 


(29) 


in the Fourier space. 


B. Intermediate friction 


For intermediate values of F, there are a considerable slip between the ambient fluid velocity v'^\z=±e and the active 
fluid velocity a = x,y. We consider the limit q T/rj' (equivalently q 2> C^), (fTSl) and (fTOl) reduce to (see 
Appendix [C| 


v'dzv'pU = -Tvp, (30) 

r)'dzVp\-e = Tvp,. (31) 


Clearly, these would be valid for wavevector qlg 1 or a system with size L Is. In addition, we should have 
pq^ <C F. Since p p'd, where d is the thickness of the 2D active fluid layer, we obtain L ^ (/gd)^^^, yielding 

Zg > L > (/gd)i/2. 

Substituting (I30p and (1311) in (fTOl) and linearising about Px = 1 and p = p^, the generalised Stokes equations for Vx 
and Vy are obtained as (see Appendix ICl) 


Vx 


V 


y 


* 2F 

* 2F 



XApqxql 
p. + < P. 

(32) 

XApqldy 

Py * 2Fg2 

(33) 

















where we have neglected rjq^Va in the limit rjq^ <C 2T. This should be valid in the wavevector range satisfying 
<C r <C q'q. With q ^ q'd, taking d 10“^m for a cortical actin layer and q' ^ N.sec/m?' for water, 
the above inequality should hold over a wide range of q. As before, Va has only active contributions at the lowest 
order in q. Similar to an ordered active polar fluid layer on a solid substrate, the hydrodynamic interactions here are 
completely cut off by the friction T and consequently v ~ 0{q) to the lowest order in the wavevector. Not surprisingly, 
Eqs. (j32l) and (|M|) are identical in structure with the form of the velocities of an active polar fluid layer resting on 
a solid surface. This is due to the fact that for qlg ^ 1, v^, i = x,y are effectively very small and hence ignored. 
This background fluid thus effectively behaves as a fixed background with the force on the 2D flow being given by 
—Fuq, (similar to a rigid substrate). Thus, with an intermediate value for T, our model active system corresponds 
surprisingly to a dry active matter, despite being in contact with an embedding bulk fluid. 

Using the above Eqs. (ElD and dSSD in eq. (1^^ . the equation for py can be written as 


dtPy 


- l)9y - (i"i + 


CA/r 1 


Dq^Py 

70 


+ iX2l2^pLqxPy - iXqyP - XxqlPy 


, ^^pqxqy / 

o ] Py n P 

r J Q . 


(34) 


Density p of the active particles still follows Eq. (l23l) . Notice that Eqs. (l34l) and (l23l) are the linearised version of the 
model Eqs. for a polar flock in Ref. [^, which is a coarse-grained model for an active polar flock in a frictional medium. 
Thus, with {IgdY/"^ <C T ^ 4, the long wavelength dynamics of our model is identical to that of a polar-ordered layer 
of a suspension of active particles on a solid substrate, an example of dry active matters. In this regime, our model 
is a representation of Ref. 0 . 


C. Weak friction limit 


In this case, T is so small that qq^ » T, or, q^ ^T/q T/{q'd) ~ l/{lsd); equivalently, L (lsd)^A_ Since d is 
small for a quasi-2D system, Ig must be very large or F very small for a physical system to display the weak friction 
limit. From the generalised Stokes Eq. for v (valid now for system size L < (ry/F)^/^) we find 


Vx = 


-i /Xppy + i^^lXppy -\- i 




Vy = -l- 


qq 

. Qqldx 

qq^ 


m 

Cq^ 

Appy + i^Appy-i 

Tfq^ Tjq^ 


Cqldx . 

-^App, 

m 

Cdldy 


App. 


(35) 

(36) 


Thus, V ~ 0{l/q) at the lowest order, in contrast to the g-dependences of the velocities for large or intermediate 
F above. The differences are due to the lack of any screening of the hydrodynamic interactions in the present case. 
Effectively, in this limit, the active fluid layer is a free standing system being completely decoupled dynamically from 
the ambient fluid. The dynamical equation for py takes the form 


X 1 

dtPy = t{— + XpAp)qyp-^ 


[Xl - l)9y - Xl + l)9x] 


AA ( 1 CApqxQy 

CAp[l-^jpy + -- P 


Dq^Py 


70 


- iX 2 ApqxPy - iXqyp - Xxq^Py 


(37) 


see, e.g., Ref. 0. Equation of motion of p is still given by Eq. 


III. LINEAR INSTABILITIES 

We now analyse the linear stability of the system from the dynamical equations obtained above by assuming a 
time-dependence for py and p of the general form exp(At). There are two independent modes, which may be static 
or moving, stable or unstable, given by two values of A. We calculate A up to the lowest order in wavevector q for 
the different cases elucidated above. 
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A. High friction limit 

Consider first strong nonequilibrium partial pressure, i.e., x/To ~ ^ ^ ApA/i. The eigenvalues A of the stability 
matrix corresponding to Eqs. (ESI) and (1^51) in polar coordinates q = {q cos 9, q sin 9), where 9 is the angle between the 
wavevector q and the ordering direction (a;-axis), up to the linear order in q are 


A = 


XX2-Kpp)An BCAfj. \ a , oaiz 

I - — -geos 6 ^ H-—geos26^ ±- {^(^2 — f^pp) cosO H-- cos20\ 

2 Sg 2 4g' 

O BKnnC 

-\-A\pno sin^ 9 + i ^ , sin 6 * sin 26* + i - 7 ^ eos6*eos20 

2 g' T]' 

-4X2KppCOS^ 6>]1/2 = A!^,A!^. 


(38) 


where, B = cos^ 9—{iyi — l) sin^ 9. Clearly, both A^,A'X scale with q and A/i. Thus, A(j:, A((: ~ q, the coefficients 

of proportionality are generally unequal and should in general be complex functions of 9 (hence anisotropic) and other 
model parameters. This linear g-dependence is different from ^-independent eigenmodes in bulk polar active fluids 
(see, e.g.. Ref. @) and is a consequence of the hydrodynamic interactions mediated by the ambient fluid. In Figs. O 
representative plots of A(ji, A((: as functions of 9 are shown for two different values of C, namely C = 1 and C = 5 for 
fixed values of other parameters and q = 1, Afi > 0. The plots clearly display a significant change in the amplitude of 
the real part with ( in one of the eigenvalues and in the amplitude of the imaginary part with C in the other one. It 
may thus be concluded that the amplitude of the unstable mode and as well as the propagating mode changes with 
change in the magnitude of C or the active stress coefficient. Figures 2] show plots of the eigenvalues as function of 9 
for negative values of C, keeping all other parameters fixed, which again shows the change in amplitude of the real and 
the imaginary parts in the two different eigenvalues, with change in value of ICj. Figures [5] compare the eigenvalues 
for different signatures of ( (the coefficient of the small fluctuations of the active stress), keeping all other parameters 
fixed. From the plots it is quite clear that the dependence of A(( and A(^: on ( is very weak. In Figs. IH the plots of 
A(( and A’X are shown for different signs and values of Ap, other parameters kept constant. The plots bring out the 
changes in the imaginary parts. Furthermore, although the real part of A(ji is always found to be positive (for these 
choices of the parameters) and signifies instability in the system, the real part of A(i shows a transition from typically 
negative values for Ap > 0 to positive values for Ap < 0. This suggests a very strong dependence of the eigenmodes 
on the value and signature of Ap or the active osmotic pressure coefficient. 

Consider in detail the limit R —0, or, 


tan^ 9 = = tan^ 9o, (39) 

ui-l 

yielding 9 = ±9o,±{9o + tt), such that B{9o) = 0. The eigenvalues in this case are 

A(0o) = i— — ^PP^^^ gcos6>o ± ^^[-(A 2 + Kpp)^cos^ 6»o -I- 4ApKo sin^ (40) 

It is evident from Fq. (IMl) that for XpKo < 0, A(0o) is fully imaginary i.e., two propagating modes, which are oppositely 
moving, are present in the system with an anisotropic q-independent wave speed. Thus, 9q gives the direction in the 
plane along which small fluctuations propagate without growth or decay at 0{q). On the other hand, for XpUg > 0 
and |4ApKo sin^ 0o| > |(A 2 + Kpp)^ cos^ 0o|, A(0o) bas a real part in addition to propagating modes. The real part 
comes from 


± [-(A 2 + Kpp)'^ cos^ 9o + AXpKo sin^ (41) 

Evidently, the real part displays instability for both signs of A/r in this case, with anisotropic decay/growth rates 
which scale with q. Thus, unlike the case with XpKo < 0, there is no particular significance of the angle 9o here. 

A schematic plot of n^Xp vs 9 for chosen values of Upp and A 2 is shown in Fig. 0 , clearly indicating the unstable 
regions and propagating modes. That the expression (14111 determines the stability at 9 = 9o be understood heuristically 
as follows. If all other parameters are set to to zero, the combination A 2 -I- Upp gives the relative speed of propagations 
of the fluctuations of Py and p in the linearised theory, with py and p being decoupled from each other. On the other 
hand the product KgAp controls the wavespeed or the growth rate of the fluctuations in the linearised coupled system 
of Py and p, depending on its sign, with all other parameters set to zero. Thus, in a situation where all the above 
four parameters are nonzero, it is generally expected that A 2 -I- Kpp has the effect of stabilising the instabilities due 
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FIG. 3: (Color online) Representative plots of the real (blue line) and imaginary (green line) parts of (a) eigenmode A(jl vs 6 
for C — 1; (b) eigenmode vs 6 for C = 1, (c) eigenmode A+ vs 6 for C = 5, and (d) eigenmode PJL vs 6 for C = 5 with fixed 
values of the other parameters; A 2 = 1, Kpp = 1/2, A/r = 1, 1/1 = 3, rj' = 1, ^ = 1, kq = 1, Ap = 2, and q = 1. Here Afj, > 0 for all 
the plots. Nonzero imaginary part implies propagating modes (see text). 


to kqXp (assuming it has the sign that corresponds to instability), by allowing reduction of local inhomogeneities to 
disperse by means of wave propagation. 

Consider the case when there are only propagating modes at 0{q) at, say, 9 = 9q. Now assume 9 very close to 9o': 
we write 9 — 9o + 69, where 69 is very small. In that case B ^ —2^i sm29o69, up to order 0{69). The eigenvalues 
corresponding to 9 = 9 q + 69 are given by 


A(0o + 69) = A{9o) 


ir]' 


g sin 26*0 cos 26*o<56* + iO{69). 


(42) 


Noting that A(0o) is fully imaginary, shows that A{9q + 69) has real parts, whose signs depend on A/r for a given 
69. Thus, we conclude that the system shows instability for either sign of A/r along with generic propagating modes 
with an anisotropic wave speed proportional to A/i. Considering A in the (g, 9) plane, we thus notice that there are 
special directions given hy 9 = ±9 q, ±(0o + t^) along which (small) perturbations move as waves without any growth 
or damping (to the linear order in g), provided A = 4ApKosin^6*o — (A 2 + Kpp)^ cos^ 9 0 < 0 is satisfied; see Fig. [S] 
Additional values of 9 for which for which the real parts of A(|i or A(/ vanish may be found from (l38l) . However, 
both the real parts will not vanish simultaneously at these angles; see Fig. [3l Along all other directions, at least one 
of A^ should have a real part, and hence perturbations will grow/decay and move. If A > 0, there are no special 
directions with only propagating modes. Since, py(x, t) and p(x, t) depend on py{q,t) and p(q, 6) for all q, hence, 
the two eigenmodes for all q, Py{:x.,t) and p(x, t) show generic moving instabilities at 0(g) for both signs of Ap for 
arbitrary choice for the active coefficients. It is also clear that at 0(g), the system can be stable only if C = 0 = C 
and XpKo < 0. Thus, the active stresses clearly distablise the system. Of course, at higher order in g, the system will 
be stabilised by large enough D or 7 pp. 
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FIG. 4: (Color online) Representative plots of the real (blue line) and imaginary (green line) parts of (a) eigenmode A(jl vs 6 
for ^ = —1, (b) eigenmode vs 6 for ^ = —1, (c) eigenmode A+ vs 6 for ^ = —5, and (d) eigenmode A(1 vs 6 for ^ = —5 
with fixed values of the other parameters; A 2 = 1, Kpp = 1/2, A/r = l,ui = 3,r]' = 1,^ = 1, kq = I, Ap = 2, and q = 1. Nonzero 
imaginary part implies propagating modes. 


It is useful to analyse the stability of the system for some particular values of 9. First we start with 0 = 0. In this 
limit the eigenvalues are given by 

A{9 = 0) = tA 2 A/rg + -iAfXKppq. (43) 

Thus there are two modes; one is purely imaginary and hence just a propagating mode, the other has both real and 
imaginary parts. The sign of the real part is determined by A/r. Thus this eigenvalue is moving and either growing 
(unstable) or decaying (stable) in time, respectively, when (i/i + l)CA/r > 0, or, < 0. 

For 0 = f, the stability eigenvalues are given by 



(z/i - l)CAfi 
877 ' 


q± 


A/rg 

2 


/ {vi - 1)C 


2 


1/2 


+ 4Ap/co 


(44) 


From Eq. dg, we note that for XpHo > 0, the system is unstable for both Afi > 0 and Afi < 0. Next, for XpKo < 0 
and (z/i — 1 )C > 0 , 

• If |4ApKo| > ^ ^ and A/i > 0, the modes are unstable and oppositely moving. 

/ _^'\ 2^2 

• However, when |4Ap«;o| < ^ ^ with A/r > 0, both the modes are unstable. There are no propagating waves. 
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FIG. 5; (Color online) Representative plots of the real (blue line) and imaginary (green line) parts of (a) eigenmode A(jl vs 6 
for C = 3i (b) eigenmode vs 6 for C = 3, (c) eigenmode A+ vs 6 for ^ = —3, (d) eigenmode A(i vs 6 for f = —3, with fixed 
values of the other parameters; A 2 = 1, Kpp = 1/2, = 1,1/1 = 3,ri' = 1,( = 3, acq = 1, Ap = 2, and q = 1. Nonzero imaginary 

part implies propagating modes. 


In the special case with </ = 0 = Ap in Eq. (1551) . i.e., if we ignore the density dependences of the active coefficients, 
the eigenvalues of the stability matrix take a simpler form 


A(C = 0 = Ap) 


iX2Afj,q COS 9 


Ari 


q cos 29, —iKppA^q cos 9, 


(45) 


which indicates the presence of propagating modes and instability for both signs of A/r above or below 9 = 

Now briefly consider the instabilities with |x/7o — A| |ApA^| (weak nonequilibrium partial pressure): Neglecting 
ApA/j, in comparison with x/lo — A, the eigenvalues A are given by 


A = \—iKpp cos9 + cos 2^ + fA2qcos6*1 ± cosff + inpp cos9]^An‘^ 

2 L ec 1 j dry' 

sin^ 9cos9 + (— — A) sin^ 9]KoAqi + AiKppAqi^ cos0[zA2 cos6* + cos20]}^/^. 


(46) 


Thus, A are no longer homogeneous functions of Ajj.. In order to progress further, assume a ’’small” Afi. Then, in an 
expansion in powers of A/i, we obtain to the lowest order in q and Afi 


A = ±,j{— — X)KoAfiqsm9, {— —A)koA/j, > 0, 


70 


7o 


- X)KoAfj,\, {—-X)koA^<0. 
7o 7o 


(47) 


A = Aiq sin 9 


(48) 
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FIG. 6: (Color online) Representative plots of the real (blue line) and imaginary (green line) parts of (a) eigenmode A+ vs 
6 for \p = 2, (b) eigenmode A(i vs 6 for Ap = 2, (c) eigenmode A+ vs 6 for Ap = 0, (d) eigenmode A(l vs 9 for \p = 0, (e) 
eigenmode A+ vs 9 for Ap = —2, (f) eigenmode A(l vs 9 for Ap = —2, with fixed values of the other parameters; A 2 = 1, ftpp = 
1/2, Ap = 1,:/! = 3, 77 ' = 1, /" = 3, Ko = 1, C = Ij q = 1. Nonzero imaginary part implies propagating modes. 


Thus, in the former case, we find instabilities for either sign of Ap, where as in the second case, we find oppositely 
moving propagating waves. 
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FIG. 7: (color online) A plot of KoXp vs Oq for some chosen values of Kpp and A 2 is shown. The regions inside the upward 
parabolas indicate the presence of moving instabilities and all other regions outside have propagating modes without damping 
or growth at 0{q) (see text). 



FIG. 8: Schematic diagram {noXp < 0 or A < 0) displaying the special angular directions given hy 6 = 9o in the plane along 
which there are only propagating modes up to 0{q) (see text). 


B. Intermediate friction 

We again consider |x/7o — A| <C |ApA^| first.To the lowest order (linear order) in q, the eigenvalues of the linear 
stability matrix are 


A = - A 2 ) cos0± l^cos^ diupp + A 2 )^ - 4koAp sin^ 9\ 

= A+,A_, (49) 

which are independent of T. We can make the following general conclusions about the mode structures from (j49p . 
First of all, none of the active stress coefficients ( and C appear in (|4^ . Thus the active stress is irrelevant in the 
dynamics to the linear order in q. The dynamics at this order in q is controlled by the remaining active coefficients, viz, 
^ppj A 2 and Ap. This is clearly in contrast to the situation with large T. Secondly, if KqXp < 0, then the discriminant 
is positive for all values of 9 to the linear order in q. Then only propagating modes will be present for all values of 9. 
As before, there should be two oppositely moving propagating modes with the speed of wave being anisotropic and 
proportional to A/r. If on the other hand KqXp > 0, the discriminant in (1491) is negative for all magnitudes of nqXp ^ 0 
at least at 0 = 7r/2, giving rise to instability in the system for both signs of A/i. These instabilities are moving in the 
opposite directions. In general, for any value of 9 satisfying cos^ 0(Kpp + A 2 )^ — d/toAp sin^ 0 > 0 and for both signs of 
A/i, both the modes are propagating without damping (or growth). Thus, any perturbation in a region of the polar 
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plane satisfying the above condition moves without any growth or decay in the amplitude (up to 0{q)). Else, in the 
remaining region of the polar plane, one of the modes is unstable and the other stable. The speed of the moving 
stabilities are unsurprisingly anisotropic and proportional to A/x. The above consideration for K^Xp > 0 allows us to 
define an angle 6 such that 

cos^ 9{Kpp + A 2 )^ — d/tgAp sin^ 0 = 0. (50) 

Then, for KgAp > 0 in the shaded region in Fig. IH] characterised by 9 there are only propagating waves at 0{q), outside 
of this region, the system is linearly unstable at 0{q). 



FIG. 9: Schematic diagram (fcoAp > 0) depicting the regions of directions in the plane (shaded region) in which only propagating 
modes exist at 0{q). Outside of the region, the system is nnstable at 0(q). In contrast, pure propagating modes without any 
damping or gworth are fonnd for Zs —t 0 (F —>■ oo) only along fonr special lines for KoAp < 0 (see text; see also Fig. [8]). 

The growth rate or relaxation rates of the unstable and stable modes are also anisotropic and scale with A/x. 
Representative plots of the real and imaginary parts of A-|_, A_ as functions of 9 for some chosen parameter values are 
shown in Fig. (nni) showing the presence of propagating modes. The regions of instabilities and propagating modes 
are clearly indicated. In particular, there are a few notable features as displayed by Fig. [101 consistent with the forms 
of the eigenvalues (H^ . For instance, for KgAp < 0, A+, A_ are wholly imaginary for all 9 and unequal, i.e., the speed 
of the two modes are different in magnitude. In contrast, for KgAp > 0, the real parts vanish over an identical range of 
9 for both the modes, that belongs to the shaded region in Fig. [HI with unequal imaginary parts, i.e., different speeds 
for the two modes. For the other values of 0, the real parts are nonzero with mutually opposite signs, representing 
stable and unstable modes, with same speeds of propagation. The overall differences with the eigenvalues for large 
(diverging) T are clearly visible. Evidently, the model is overall stable at the linear order in q, provided, ApKg < 0- 
In this stable sector of the parameter space, the results of Ref. Q that includes the effects of the nonlinearities and 
noises should directly apply here. Lastly, dry active matters are characterised by density segregation in the steady 
states [11,14^. Our linearised treatment is unable to capture this. 

We now consider briefly the case with lx /70 — A| |ApA^|. Proceeding as in Sec. IIII Al above, the eigenvalues to 
the lowest order in q and A/x are given by 

A = ±,j {—— A)KoA/xgsin0, (51) 

V 7o 

yielding instability for (x/70 ~ A)ko A/x > 0 and oppositely moving propagating modes for (x/70 “ A)ko A/x < 0. These 
results are identical to the corresponding results in Sec. IIII Al 


C. Weak friction limit 


We now analyse the linear instabilities for rjq^ ^ T. One of the eigenvalues A of the linear stability matrix is 
non-zero at 0{q^). We find 


A = 


CA/i 

20 


(cos^ 0 — sin^ 0 ) \{vi — 1 ) sin^ 0 — (x^i -|- 1 ) cos^ 


(52) 
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FIG. 10: (Color online) Representative plots of the real (thin black line) and imaginary (thick line) parts of A+ and A_ 
(calculated up to 0{q)) vs 6 for some chosen parameter values A 2 = 1, Kpp = 2, Afi = 1 and q — 1. (a) A^. vs 6 for KoXp = 4, 
(b) A- vs 6 for KoXp = 4, (c) A+ vs 0 for KoXp = —4, and (d) A_ vs 0 for KoXp = —4. In (a) the real part (thin line) is located in 
the region of negative A+ indicating stability and the imaginary part (thick line) indicates the presence of propagating modes. 
In (c) and (d) there are no real parts (black thin line) when KoXp < 0. Unlike the case for U —>■ 0, there are ranges of directions 
in the plane where there are only propagating modes at 0(q) (see text). 


With a given choice for the sign of CA/r (say positive), A > (<)0 for cos^ 0 — sin^ 0 and (z/i — 1) sin^ 0 — (i^i + I) cos^ 0 
having the same (opposite) signs and vice versa for (Afi < 0, suggesting instabilities for either sign of (Afx. These 
results are identical to those in Ref. @ for a bulk polar ordered active fluid. It is not a surprise that our results are 
same as those in Ref. for in the weak T limit, the active fluid layer in our model is effectively dynamically decoupled 
from the ambient fluid and hence acts as a free standing system, and hence, identical to the system considered in 

Ref. 0. 


IV. LINEAR INSTABILITIES AND MEASUREMENTS OF F 

As our results above reveal, the magnitude of F delineates different regimes of the model. While all these regimes 
display generic long wavelength instabilities in the different regions of the parameter space, the detailed nature of 
the instabilities and the regions in the parameter space where they are present, vary depending on F. For easy 
comparison, we provide here a table (Table H]) which differentiates between the instabilities in the three different 
regimes, as delineated by T (assume lx /70 — A| <C |ApA^|): 

Despite the loose similarities between the nature of the long wavelength instabilities for large and moderate T, closer 
inspection reveals significant differences between the two cases. With a large (formally diverging) T (Is/L <C 1), the 
active fluid velocities Va ~ 0(q°), where as, for moderate T, Va ~ 0(q), a = x,y. Furthermore, with a diverging F, 
the system is unstable in the full parameter space to the lowest order in q along all angles in the polar plane, except 
for along the lines 0 = ±0O)±(^o + t’")- Along these special directions, there are only propagating modes without 
any damping or growth (to the linear order in q). At every other value of 0, one mode is unstable. Thus, for all 
(finite) values of the parameters and both signs of A/r, there are moving instabilities with anisotropic speeds. In 
contrast, with intermediate interfacial friction {Ig/L ^ 1), there are regions in the parameter space where there are 
only propagating modes with no instabilities at the lowest order in q for any 0; only in a subspace of the parameter 
space, one encounters moving instabilities for either sign of Ay. Even in such a parameter subspace, there are only 
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High friction (L ^ Zg) 

Intermediate friction (L <C Zs) 

Weak friction [L (Zgd)^/^] 

Eigenvalues vanish at 0(g°). 

Again eigenvalues vanish at O(q^). 

Nonzero eigenvalues: generic linear 
instabilities for at 0{q°) for both 
signs of A/r; no propagating modes. 

0{q): Generically linearly unstable at 0{q) for 
both signs of Aq. In the plane at an angle 9o 
(measured with respect to the direction of the 
reference orientation) given by tan^ 9o = {vi + 
l)/(i/i — 1) for ApAco < 0, or, for ApAco > 0 and 
4ApAvosin^0o < (A 2 + Atpp)^ cos^ 0o: only propa¬ 
gating modes at 0{q). Else instability even at 
9 = 9o. 

0(q): Only propagating modes for AtoAp < 
0. Else, for AtoAp > 0, there are only 
propagating modes without any damping 
or growth for angles 9 in the plane satis¬ 
fying cop 9 {Kpp A 2 )^ — 4KoApSin^6' > 0 
and instability elsewhere in the plane. 



TABLE I: Table describing the instabilities in different frictional regimes. 


propagating modes at 0{q) for a range of 9; for other values of 9, moving instabilities are present. 

At this stage, it is useful to compare with available experimental results. To do this, numerical estimates of the slip 
length Is or the slip coefficient T are needed. To our knowledge, systematic measurements of T or Ig for active fluids 
are lacking. Nonetheless, based on the available information we can make the following comments. Ref. reports 
a hydrodynamic length I = (ry/T)^/^ to be of the order of 10 /rm. Assuming rj ^ r/'d and d 200nm as the 
thickness of actin cortex, we find Ig ^ P/d ^ lO^/im, larger than the typical size of cortical action layers (~ 100/rm). 
Nevertheless, our estimation of Ig is not precise and hence it is difficult to comment upon the experimental realisability 
of the high friction case of in-vitro cortical actin layers in water based upon our estimates. On the other hand, L 
smaller than Ig should correspond to the intermediate friction case; we expect this to be realised in experiments on 
cortical actin layers in water. Lastly, for L < {Igd^P ^ 10~^fim, which is certainly small, the system should behave 
as a free standing system. Equivalently, for a larger system with L > the dynamics of a free standing 

film will be observed if the system is probed at length scales much smaller than ~ see, e.g.. Ref. [^ . 

The sensitive dependences of the long wavelength dynamics of the model on T may be used to make experimental 
estimates of T in a given system (i.e., for a fixed values of all other parameters including L). We note that directly 
connecting our theoretical predictions with experimental results is not an easy task. Nonetheless, given the generic 
nature of our continuum active fluid theories and noting that since an experiment is necessarily performed on finite 
systems, long wavelength limit should imply 27r/g —>■ L, where g is a wavevector of interest and L the system size, 
it is expected that all experiments that may be described by the same long wavelength continuum equations should 
display similar long wavelength linear instabilities, characterised by their growth rates or thresholds of the linear 
instabilities. In general, these properties should allow us to compare the theoretical predictions with experiments, 
at least qualitatively, although varying L experimentally is expected to be a challenging task. More specifically, we 
can make the following comments. First of all, if linear instabilities are found to persist for all 9, then our above 
results indicate that Ig/L ^ 1 or a large T: T d'/h, hence, T setting a lower bound for T ~ r]'/L. 

On the other hand, if linear instabilities are not found at 0{q) or found only over a range of 9, we can conclude 
Ig/L ^ 1, or, r <C Tq'/L^ giving an upper bound on T. At the same time, we must have rjq^ « T,ri ^ rj'd, yielding 
r rj'd/a'^, where a is a small scale molecular cut off), such that for a~^ < q, the continuum theory breaks 
down. This provides a lower limit on T. In contrast, the system behaves as a free standing 2D film for L <C {rj/VyP. 
Information on T may also be obtained by measuring the correlation functions of the local velocity fields, orientation 
and density fluctuations and using the relations p2p and (1331) when the system is stable with an intermediate T, i.e., 
with KoXp < 0. Velocity helds may be measured, e.g., by attaching a small bead with the actin hlament and tracking 
its instantaneous positions. Orientation and density fluctuations may be measured by optical methods and scattering 
experiments, respectively. Lastly, as we have shown in Appendices (ID II) . (|D 21) and (ID 31) that the magnitude of bulk 
3D fluid velocity depends strongly on T for q'/L T q'd/a?. Thus, measurement of the ambient 3D velocity field, 
e.g., by tracking the position of a tracer particle, should also be helpful in extracting numerical estimates on T. 


V. NEMATIC LIMIT OF THE DYNAMICS 

Until now we have considered polar active particles, so that the corresponding dynamics is not invariant under 
p —^ — p. In the nematic limit, the dynamics is invariant under p —>■ — p. Hence, active coefficients X 2 ,Xp,Ko and Kpp 
and equilibrium couplings y and A are zero. Thus, to the lowest order, the dynamical equations in the strong friction 
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case are 


dpy 




dt Arj'q 

dp 




CApq^qy 
Py + - -o - P 


^ = -Jppq'^P + wApqlp + wpoApq^qyPy. 


As before, assume a time dependence of the form exp(At) for the fluctuations. Then, to the lowest order in q 

A _ 0, Q _ gjj^2 0^ _ 2^^ gjj^2 0 _|_ (,pg2 _ 

4:7]' 


(53) 

(54) 

(55) 


Thus, with positive A is positive (negative) for cos^ 9 — sin^ 6 and — 1) sin^ 9—{i'i + 1) cos^ 9 having the same 
(opposite) signs. Similarly for CA/r < 0. Thus, the system is unstable for both signs of Ap. 

For a finite T, to the lowest order in q, the corresponding dynamical equations with nematic symmetry are 


~'~Py ~ Py^ 

-IppQ^P + wApqlp + wpoApq^qyPy. (57) 

Interestingly, Eqs . (IMl) and (EH) are identical to those in Ref. M for active nematics on a substrate. Thus, the 
results of Ref. [1^ are to hold here. We do not discuss these here in details. Regardless of the details, in the nematic 
limit there are no propagating waves and the instabilities are always static or localised. In contrast, active polar 
ordered systems are characterised by the presence of generic propagating modes and moving instabilities. Finally, the 
eigenmodes in both the nematic and polar ordered systems with strong interfacial friction with the embedding fluid 
scale with q. However, for intermediate friction, the eigenmodes for the nematic system scale as q^, where as for the 
corresponding polar ordered system, they scale as q. 


dPy 

dt 

dt 


VI. SUMMARY 

In this work, we have set up the generic coarse-grained dynamics of a thin layer of polar ordered active particle 
suspensions frictionally coupled to the bulk isotropic passive fluid with an arbitrary friction coefficient F. In a linearised 
treatment for small fluctuations around uniformly polar ordered states, we show that our model describes a layer of 
wet active matter, dry active matter and a free standing film, respectively, for L ^ p'/T, r/'/F ^ L ]S> {p'd/ry^"^ and 
p'd/r ^ L^. The nature and the conditions for linear instabilities in the long wavelength limit depend sensitively on 
F. These features may be used to find estimates about F in a given 2D active fluid layer embedded in a bulk passive 
fluid. We also discuss the nematic limit of the dynamics and compare it with their polar analogues. 

Our results evidently highlight the crucial role played by the interfacial friction and demonstrate how experimental 
knowledge about the linear instabilities may be used to extract information about the friction coefficient. Actual 
biological realisations of quasi-2D active fluids have more complicated structures. Our work should be considered 
only as a first step towards a more complete physical understanding of such systems. We expect our results to be 
useful in understanding in-vitro experiments on reconstituted layers of ordered actin filaments with molecular motors 
in an embedding fluid (e.g., water). Experimental validation of the F-dependences of the linear instabilities are 
expected to be highly challenging tasks. Nevertheless, we look forward to possible experimental attempts to study 
the issues highlighted here. Lastly the formal similarities between the dynamical equations with moderate interfacial 
friction and those for a polar ordered system resting on a solid substrate open up the possibilities of studying the 
physics of moderate friction by performing experiments on an analogous system resting on a solid substrate. 

Our analyses are valid for small fluctuations around an ordered state. Thus no conclusions may be drawn from 
our studies about the eventual steady states in the event of the linearly unstable uniform initial states. Numerical 
solutions of the full model equations should yield valuable information in this regard. We made several simplifying 
assumptions while setting up our framework. For instance, we have assumed the active fluid layer to be inflexible and 
hence the out-of-plane fluctuations are prohibited. However, this condition may be violated for reconstituted actin 
filaments on a liposome. Thus for better quantitative understanding of the experimental results, a thin layer of active 
fluid with finite flexibility (i.e., with a finite surface tension or bending modulus) should be studied. Secondly, the 
system may not even be overall flat and may have a finite curvature. In this case, our results should hold over scales 
smaller than the radius of curvature. Our assumption of equal friction coefficients on both the sides of the active 
system is also a simplification. Generalisation to unequal frictions on both sides may be done in a straightforward 
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way. It will be interesting to stndy the diffusivity of a test particle inside the 2D active polar system. It is well-known 
that the diffusivity of a test particle in a free standing thin active fluid layer shows starkly unusual properties, e.g ., 
dependences on the thickness [i^, in contrast to the diffusivity of a small particle in a quasi-2D passive fluid [d^. 
Given our results here, it is expected that the diffusivity of a test particle in a 2D polar ordered medium is affected 
by the interplay of hydrodynamic interactions by the embedding medium and the strength of the interfacial friction. 
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Appendix A: Derivation of the full 2D generalised Stokes equation for Vi 


Here, we derive the full 2D generalised Stokes equation. Using Eqs. ([7]), (ITUl) in Eq. (O, the generalised Stokes Eq. 
may be written as 

+ Aflda{C{p)PaPp) + '^daiPahp + ppha) - ^^aiPadfiPp + PfjdaJIp) + ^da{Pahp - ppha) = dpW - F>(A1) 
Using incompressibility H can be derived from m as 

n = A/r^^(C'(p)(p«p^) -f y^^(pa^/3 +P;3/la) 

-^^^{Apad/sp + AppdaP + XPadpV ■ p -f XPpdaV ■ p) -f (A2) 

where 1/V^ is the inverse of V^. 

Using this value of H in (jAll) . the Stoke’s equation is derived as 






ApPp^doL\(^ (p)Pq,P^] J- ^ Ppjda{pahj J- Pjhoi') ^oPpjda{Ap(^d^p -t- Ap^df^p -t- xPa^j^ ■ P T XP^^ol^ ■ p) 

T“F/57^q (Pa^7 P’y^Oi) — (-^^) 


where Pap is the transverse projection operator written as Pap = Sap — . Let us redefine PjpFj as Fp given by 

Linearising about = 1, the Stokes eqnation (|A3I) is simplified to 


pV Vj^ -\- Pa\ [C P7] 2 2 

-A^Pp^d^d^p - x^Pp^d.d^V ■ P - A^Pp^dlp - XjPp.dlV ■ P = -Fp, (A4) 

Now using Eq. m. we find hy = — = DV^py + xdyP- In addition, hx acts as a Lagrange multiplier to enforce the 

constraint p^ = 1. Notice that hy contributes terms which are higher order in gradients in Eq. (TAd)) . Thus neglecting 
all the higher order terms, the generalised Stokes equation up to the lowest order in gradients is given by Eq. dm). 


Appendix B: Fp for high friction {L 3> L) 

The velocity and hydrodynamic pressure for the subphase and superphase are given by Eqs. (IE5])-(IBT^. We impose 
incompressibility on the 3D ambient fluid: 


dW. = - 


(Bl) 


Viu' with i=x,y. 
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for both z > 0 and z < 0. Fourier transforming the in-plane coordinates x = {x,y), 

+ = 9.n', 

r]'{—q‘^ + dl)vl = iqiW, and 

(_g 2 ^ ^ 2 )^' = 0 , 

where z = a;, y; q = {qx, qy) is the in-plane Fourier wavevector. The above equations can be solved 
the solutions for Vy, v'^ and 11'. We write 

v'x = {Ai + Biz) exp(-gz) for z > 0 , 

= {A 2 + B 2 Z) exp {qz) for z < 0 , 
v'y = (A 3 -I- B 3 Z) exp {-qz) for z > 0, 

= (A 4 -I- B^z) exp {qz) for z < 0, 
v'^ = (Cl -I- I?iz) exp {—qz) for z > 0, 

= (C 2 -I- D 2 Z) exp {qz) for z < 0, and 
n' = El exp {—qz) for z > 0 , 

= E 2 exp {qz) for z < 0 , 

where coefficients Ai, A 2 ,C 2 are real or imaginary functions of q. 

The incompressibility condition m yields 


(B2) 

(B3) 

(B4) 

together to obtain 

(B5) 

(B6) 

(B7) 

(B8) 

(B9) 

(BIO) 

(Bll) 

(B12) 


Di = -iqxAi - iqyAs + qCi = i — Bi + i—Bs, 

q q 

D 2 = -iqxA2 - iqyAi - qC2 =-i—B2 - i — B^. 

q q 


The continuity of velocity or Eq. (HI gives 


(B13) 
(B14) 

(B15) 
(B16) 
(B17) 

As the active fluid film is two dimensional, there is no discontinuity over the vertical gradient of v'^ (since Vz =0). 
This allows us to write 


Ai = A2 = Vx 

A3 = A4 = Vy 

Cl = C 2 


^zXz\z—e — ^zX z\z— — e-i 


(B18) 


which yields using Eqs. (IB9|) and (jBlOl) 

Di=2qCi+D2. (B19) 

The tangential stress Fi may be evaluated using the Stokes equation dm. The Stokes equation for v{ yields 

ry'V^V3 X u' = 0. (B20) 

Eq. (IB20I) gives us further relations 

(B21) 


B 3 = -i^Di = ^Bi and 


Qx 

B4 = i^D2 = ^B2. 
q qx 

Using Eqs. (IBT^ - (IBT711 . dHH), (IB2T]) and (IB22I) . the x-component of 3D force F is obtained as 
Fx = r]'{dzv'x + dxv'Jle - r]'{dzv'x + dxv'^)\-e 


(B22) 


'2 

= r]'[—2—Vx — 2.qVx — 2^^—^Vy -\- i — {Di + D 2 ) + Bi — i? 2 ] 


q 


= r]'[- 2 —Vx- 2 qvx- 2 ^^^Vy\ 


= -2ri'qvx, (B23) 

where we have used incompressibility of the 3D velocity in the last line. Similarly we get the y-component of F as 

Fy = -2r]'qvy, (B24) 

using the no-slip condition equating v with the in-plane components of v' at the active fluid layer. 
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Appendix C: Form of Fp for intermediate friction {h ^ L ^ 


We start with 


dv' 

V'-^\z=e = T^iVa\z=e -Va), a = X,y. 


(Cl) 


A similar condition exists at z = — e. Now using the forms of v'^, v' and as given by (lB5l) . (lB7l) and (IB9I) we obtain 


r]'{-Aiq + Bi) = r(Ai - v^), 

+ Bs) = r(A 3 - Vy). 

In the weak friction limit, F ^ 0{ri'q) in the wavevector range of interest. Thus, 

r]'{-Aiq + Bi) = 

r]'{-A3q + B3) =T{-Vy). 


or, equivalently. 


' - _r 

V — FUa;, 

dv' 


at z = e. Similar considerations at z = —e finally yields 




dz 


This yields for the 2D generalised Stokes equation which Va satisfy 

TyV Uq, + T C^pBax^yPy T C^pBax^xP ‘2YVa — 0, . 

Now write Eq. (IC9I) in the Fourier space and neglect pq^Va assuming -qq^ <§; 2F. This yields Eqs. ([32]) and dMl). 


(C2) 

(C3) 


(C4) 

(C5) 


(C6) 

(C7) 

(C8) 

(C9) 


Appendix D: Velocity profiles of the ambient fluid 

It is instructive to obtain the flow profiles of three-dimensional velocity fields, that are created by the (small) 
fluctuations in p and Py, in the three different regimes of our model as delineated by the values of F. 


1. Large F (L » 1 ^) 

In this case v^{x,y,z = ±e) = Vi{x,y), i = x,y. Since n(,(z = ±e) = 0, from (jB9l) and (IBIOI) . Ci = 0 = € 2 - Using 
the no-slip condition on n'(z = ±) and the 3D incompressibility oi v'^,a = x, y, z, {iqxv'x + *9y)|z=±e = 0 = ^^\z=± in 
the Fourier space. This yields Di = 0 = D 2 . Thus, = 0 everywhere above and below the active fluid layer. Hence, 
the flow in the surrounding fluid is actually 2D, parallel to the active fluid layer. We further find Bi = 0 = B 2 and 
B 3 = 0 = B 4 . Thus in the Fourier space, 

v'i{qx,qy,z) = v,{qx,qy)exp{-qz), z > 0, 

= v^{qx,qy)exp{qz), z < 0. (Dl) 

Therefore, n' has the same form as Vi with an exponentially damped amplitude by a factor exp(—g|z|) and hence 
shows the same instabilities at 0 {q). 
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2. Intermediate F (Zs S> i 3> (Isd)^^^) 


In the intermediate friction case, the 3D shear stress balance is given by 

^dv', 


Using the above equations and 


V 


,H\ 

dz '■ 


= -bUile, 


(D2) 

(D3) 


and (IB8I) . we get a set of relations between the couplings given by 


-qAi+Bi = - -Va:, 

1 ] 

r 

qA2 + B 2 = —Vx, 

V 

T 

-qA3 + B3 = - -Vy, 


qA^ + Bi = —Vy. 

V 


(D4) 

(D5) 

(D6) 

(D7) 


Now using (lD4l) . (IDSI) . (IDel) . (iDTl) . 3D incompressibility (V-v' = 0) of ambient fluid, 2D incompressibility (Vj_ ■ v = 
0) of active fluid layer and equations (IB21|) and (IB22|) . we can show that 


Bi = B 2 = B 3 = Bi = 0 and, 
Di = D 2 = 0 . 


(D8) 

(D9) 


Hence we obtain, = 0 for all z identically and v'j^{qx,qy, z) = Tvi exp{—q\z\)/{rj'q) and has the same dependences 
on p and py as for large F. Thus, v'^, a = x,y, z is again 2D. Nonetheless, u' for moderate F is different from u' for 
large F, since the solutions for p and Py have very different explicit forms for moderate F. In particular, v[{x^y,z) 
shows instability only if K^Xp > 0, the same condition for instability for the 2D active fluid layer. Notice that in the 
present case in addition to the exp(—qlzj) factors, v[ is further scaled down in comparison with Vi by a factor F/ (rj'q), 
and hence should be small in magnitude. 


3. Small r (L^ < Isd) 

In the limit of very small F, we have ri'^^\z=± ~ 0. i.e., we effectively have the stress-free boundary condition on 
u■ at 2 : = 0. In addition, = 0 at z = 0. Now using the results in Sec. O we find Ai = A 2 = Bi = B 2 = A 3 = A 4 = 
B 3 = B 4 = 0, i.e., z;' = 0 at all z > 0 and z < 0 identically. In addition, = 0 everywhere. Thus, the 3D velocity 
field vanishes. 
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